TAD boundary deletion causes PITX2-related cardiac electrical and structural defects.

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity.

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

RPL35A drives ovarian cancer progression by promoting the binding of YY1 to CTCF promoter.

Ovarian cancer is one of the most common gynaecological malignancies with poor prognosis and lack of effective treatment. The improvement of the situation of ovarian cancer urgently requires the exploration of its molecular mechanism to develop more effective molecular targeted drugs. In this study, the role of human ribosomal protein l35a (RPL35A) in ovarian cancer was explored in vitro and in vivo. Our data identified that RPL35A expression was abnormally elevated in ovarian cancer. Clinically, high expression of RPL35A predicted short survival and poor TNM staging in patients with ovarian cancer. Functionally, RPL35A knock down inhibited ovarian cancer cell proliferation and migration, enhanced apoptosis, while overexpression had the opposite effect. Mechanically, RPL35A promoted the direct binding of transcription factor YY1 to CTCF in ovarian cancer cells. Consistently, RPL35A regulated ovarian cancer progression depending on CTCF in vitro and in vivo. Furthermore, RPL35A affected the proliferation and apoptosis of ovarian cancer cells through PPAR signalling pathway. In conclusion, RPL35A drove ovarian cancer progression by promoting the binding of YY1 and CTCF promoter, and inhibiting this process may be an effective strategy for targeted therapy of this disease.

Dissection of a CTCF topological boundary uncovers principles of enhancer-oncogene regulation.

Enhancer-gene communication is dependent on topologically associating domains (TADs) and boundaries enforced by the CCCTC-binding factor (CTCF) insulator, but the underlying structures and mechanisms remain controversial. Here, we investigate a boundary that typically insulates fibroblast growth factor (FGF) oncogenes but is disrupted by DNA hypermethylation in gastrointestinal stromal tumors (GISTs). The boundary contains an array of CTCF sites that enforce adjacent TADs, one containing FGF genes and the other containing ANO1 and its putative enhancers, which are specifically active in GIST and its likely cell of origin. We show that coordinate disruption of four CTCF motifs in the boundary fuses the adjacent TADs, allows the ANO1 enhancer to contact FGF3, and causes its robust induction. High-resolution micro-C maps reveal specific contact between transcription initiation sites in the ANO1 enhancer and FGF3 promoter that quantitatively scales with FGF3 induction such that modest changes in contact frequency result in strong changes in expression, consistent with a causal relationship.

Tet-mediated DNA methylation dynamics affect chromosome organization.

DNA Methylation is a significant epigenetic modification that can modulate chromosome states, but its role in orchestrating chromosome organization has not been well elucidated. Here we systematically assessed the effects of DNA Methylation on chromosome organization with a multi-omics strategy to capture DNA Methylation and high-order chromosome interaction simultaneously on mouse embryonic stem cells with DNA methylation dioxygenase Tet triple knock-out (Tet-TKO). Globally, upon Tet-TKO, we observed weakened compartmentalization, corresponding to decreased methylation differences between CpG island (CGI) rich and poor domains. Tet-TKO could also induce hypermethylation for the CTCF binding peaks in TAD boundaries and chromatin loop anchors. Accordingly, CTCF peak generally weakened upon Tet-TKO, which results in weakened TAD structure and depletion of long-range chromatin loops. Genes that lost enhancer-promoter looping upon Tet-TKO showed DNA hypermethylation in their gene bodies, which may compensate for the disruption of gene expression. We also observed distinct effects of Tet1 and Tet2 on chromatin organization and increased DNA methylation correlation on spatially interacted fragments upon Tet inactivation. Our work showed the broad effects of Tet inactivation and DNA methylation dynamics on chromosome organization.

DARDN: A Deep-Learning Approach for CTCF Binding Sequence Classification and Oncogenic Regulatory Feature Discovery.

Characterization of gene regulatory mechanisms in cancer is a key task in cancer genomics. CCCTC-binding factor (CTCF), a DNA binding protein, exhibits specific binding patterns in the genome of cancer cells and has a non-canonical function to facilitate oncogenic transcription programs by cooperating with transcription factors bound at flanking distal regions. Identification of DNA sequence features from a broad genomic region that distinguish cancer-specific CTCF binding sites from regular CTCF binding sites can help find oncogenic transcription factors in a cancer type. However, the presence of long DNA sequences without localization information makes it difficult to perform conventional motif analysis. Here, we present DNAResDualNet (DARDN), a computational method that utilizes convolutional neural networks (CNNs) for predicting cancer-specific CTCF binding sites from long DNA sequences and employs DeepLIFT, a method for interpretability of deep learning models that explains the model's output in terms of the contributions of its input features. The method is used for identifying DNA sequence features associated with cancer-specific CTCF binding. Evaluation on DNA sequences associated with CTCF binding sites in T-cell acute lymphoblastic leukemia (T-ALL) and other cancer types demonstrates DARDN's ability in classifying DNA sequences surrounding cancer-specific CTCF binding from control constitutive CTCF binding and identifying sequence motifs for transcription factors potentially active in each specific cancer type. We identify potential oncogenic transcription factors in T-ALL, acute myeloid leukemia (AML), breast cancer (BRCA), colorectal cancer (CRC), lung adenocarcinoma (LUAD), and prostate cancer (PRAD). Our work demonstrates the power of advanced machine learning and feature discovery approach in finding biologically meaningful information from complex high-throughput sequencing data.

STAG2: Computational Analysis of Missense Variants Involved in Disease.

The human STAG2 protein is an essential component of the cohesin complex involved in cellular processes of gene expression, DNA repair, and genomic integrity. Somatic mutations in the STAG2 sequence have been associated with various types of cancer, while congenital variants have been linked to developmental disorders such as Mullegama-Klein-Martinez syndrome, X-linked holoprosencephaly-13, and Cornelia de Lange syndrome. In the cohesin complex, the direct interaction of STAG2 with DNA and with NIPBL, RAD21, and CTCF proteins has been described. The function of STAG2 within the complex is still unknown, but it is related to its DNA binding capacity and is modulated by its binding to the other three proteins. Every missense variant described for STAG2 is located in regions involved in one of these interactions. In the present work, we model the structure of 12 missense variants described for STAG2, as well as two other variants of NIPBl and two of RAD21 located at STAG2 interaction zone, and then analyze their behavior through molecular dynamic simulations, comparing them with the same simulation of the wild-type protein. This will allow the effects of variants to be rationalized at the atomic level and provide clues as to how STAG2 functions in the cohesin complex.

Differential DNA methylation and CTCF binding between the ESR1 promoter a of MCF-7 and MDA-MB-231 breast cancer cells.

BACKGROUND: ESR1 is expressed by 60-70% of breast tumours. it's a good prognosis factor and the target of hormone therapy. Optimization of ESR1 reactivation therapy is currently ongoing. Here we probe if the transcription factor CTCF plays a role in the differential expression of ESR1 in the breast cancer cell lines MCF-7 (ESR1+) and MDA-MB-231 (ESR1-). METHODS AND RESULTS: Knockdown of CTCF in MCF-7 resulted in decreased ESR1 gene expression. CTCF binds to the promoter of ESR1 in MCF-7 but not in MDA-MB-231 cells. CTCF ESR1 binding sites are unmethylated in MCF7 but methylated in MDA-MB-231 cells. CONCLUSION: ESR1 expression in MCF7 cells is dependent on CTCF expression. CTCF can bind to specific regions of the promotor of ESR1 gene in MCF-7 cells but not in MDA-MB-231 cells, this correlates with the methylation status of these regions and could be involved in the transcriptional regulation of ESR1.

CTCF regulates hepatitis B virus cccDNA chromatin topology.

Hepatitis B Virus (HBV) is a small DNA virus that replicates via an episomal covalently closed circular DNA (cccDNA) that serves as the transcriptional template for viral mRNAs. The host protein, CCCTC-binding factor (CTCF), is a key regulator of cellular transcription by maintaining epigenetic boundaries, nucleosome phasing, stabilisation of long-range chromatin loops and directing alternative exon splicing. We previously reported that CTCF binds two conserved motifs within Enhancer I of the HBV genome and represses viral transcription, however, the underlying mechanisms were not identified. We show that CTCF depletion in cells harbouring cccDNA-like HBV molecules and in de novo infected cells resulted in an increase in spliced transcripts, which was most notable in the abundant SP1 spliced transcript. In contrast, depletion of CTCF in cell lines with integrated HBV DNA had no effect on the abundance of viral transcripts and in line with this observation there was limited evidence for CTCF binding to viral integrants, suggesting that CTCF-regulation of HBV transcription is specific to episomal cccDNA. Analysis of HBV chromatin topology by Assay for Transposase Accessible Chromatin Sequencing (ATAC-Seq) revealed an accessible region spanning Enhancers I and II and the basal core promoter (BCP). Mutating the CTCF binding sites within Enhancer I resulted in a dramatic rearrangement of chromatin accessibility where the open chromatin region was no longer detected, indicating loss of the phased nucleosome up- and down-stream of the HBV enhancer/BCP. These data demonstrate that CTCF functions to regulate HBV chromatin conformation and nucleosomal positioning in episomal maintained cccDNA, which has important consequences for HBV transcription regulation.

Transcriptional regulation of FACT involves Coordination of chromatin accessibility and CTCF binding.

Histone chaperone FACT (facilitates chromatin transcription) is well known to promote chromatin recovery during transcription. However, the mechanism how FACT regulates genome-wide chromatin accessibility and transcription factor binding has not been fully elucidated. Through loss-of-function studies, we show here that FACT component Ssrp1 is required for DNA replication and DNA damage repair and is also essential for progression of cell phase transition and cell proliferation in mouse embryonic fibroblast cells. On the molecular level, absence of the Ssrp1 leads to increased chromatin accessibility, enhanced CTCF binding, and a remarkable change in dynamic range of gene expression. Our study thus unequivocally uncovers a unique mechanism by which FACT complex regulates transcription by coordinating genome-wide chromatin accessibility and CTCF binding.

SENP1 knockdown-mediated CTCF SUMOylation enhanced its stability and alleviated lipopolysaccharide-evoked inflammatory injury in human lung fibroblasts via regulation of FOXA2 transcription.

BACKGROUND: Excessive inflammation is the main cause of treatment failure in neonatal pneumonia (NP). CCCTC-binding factor (CTCF) represents an important node in various inflammatory diseases. In the present study, we tried to clarify the function and underlying molecular mechanism of CTCF on an in vitro cellular model of NP, which was generated by simulating the human lung fibroblast cell line WI-38 with lipopolysaccharide (LPS). METHODS: The SUMOylation level and protein interaction were verified by Co-immunoprecipitation assay. Cell viability was measured by Cell Counting Kit-8 assay. Inflammatory factors were examined by Enzyme-linked immunosorbent assay. Cell apoptosis was evaluated by TUNEL assay. The binding activity of CTCF to target promoter was tested by chromatin immunoprecipitation and luciferase reporter assay. RESULTS: LPS treatment restrained cell viability, promoted the production of inflammatory factors, and enhanced cell apoptosis. CTCF overexpression played anti-inflammatory and anti-apoptotic roles. Furthermore, CTCF was modified by SUMOylation with small ubiquitin-like modifier protein 1 (SUMO1). Interfering with sumo-specific protease 1 (SENP1) facilitated CTCF SUMOylation and protein stability, thus suppressing LPS-evoked inflammatory and apoptotic injuries. Moreover, CTCF could bind to the forkhead box protein A2 (FOXA2) promoter region to promote FOXA2 expression. The anti-inflammatory and anti-apoptotic roles of CTCF are associated with FOXA2 activation. In addition, SENP1 knockdown increased FOXA2 expression by enhancing the abundance and binding ability of CTCF. CONCLUSIONS: SUMOylation of CTCF by SENP1 knockdown enhanced its protein stability and binding ability and it further alleviated LPS-evoked inflammatory injury in human lung fibroblasts by positively regulating FOXA2 transcription.

Outward-oriented sites within clustered CTCF boundaries are key for intra-TAD chromatin interactions and gene regulation.

CTCF plays an important role in 3D genome organization by adjusting the strength of chromatin insulation at TAD boundaries, where clustered CBS (CTCF-binding site) elements are often arranged in a tandem array with a complex divergent or convergent orientation. Here, using Pcdh and HOXD loci as a paradigm, we look into the clustered CTCF TAD boundaries and find that, counterintuitively, outward-oriented CBS elements are crucial for inward enhancer-promoter interactions as well as for gene regulation. Specifically, by combinatorial deletions of a series of putative enhancer elements in mice in vivo or CBS elements in cultured cells in vitro, in conjunction with chromosome conformation capture and RNA-seq analyses, we show that deletions of outward-oriented CBS elements weaken the strength of long-distance intra-TAD promoter-enhancer interactions and enhancer activation of target genes. Our data highlight the crucial role of outward-oriented CBS elements within the clustered CTCF TAD boundaries in developmental gene regulation and have interesting implications on the organization principles of clustered CTCF sites within TAD boundaries.

Transcription induces context-dependent remodeling of chromatin architecture during differentiation.

Metazoan chromosomes are organized into discrete spatial domains (TADs), believed to contribute to the regulation of transcriptional programs. Despite extensive correlation between domain organization and gene activity, a direct mechanistic link is unclear, with perturbation studies often showing little effect. To follow chromatin architecture changes during development, we used Capture Hi-C to interrogate the domains around key differentially expressed genes during mouse thymocyte maturation, uncovering specific remodeling events. Notably, one TAD boundary was broadened to accommodate RNA polymerase elongation past the border, and subdomains were formed around some activated genes without changes in CTCF binding. The ectopic induction of some genes was sufficient to recapitulate domain formation in embryonic stem cells, providing strong evidence that transcription can directly remodel chromatin structure. These results suggest that transcriptional processes drive complex chromosome folding patterns that can be important in certain genomic contexts.

Antigen exposure reshapes chromatin architecture in central memory CD8+ T cells and imprints enhanced recall capacity.

CD62L+ central memory CD8+ T (TCM) cells provide enhanced protection than naive cells; however, the underlying mechanism, especially the contribution of higher-order genomic organization, remains unclear. Systematic Hi-C analyses reveal that antigen-experienced CD8+ T cells undergo extensive rewiring of chromatin interactions (ChrInt), with TCM cells harboring specific interaction hubs compared with naive CD8+ T cells, as observed at cytotoxic effector genes such as Ifng and Tbx21. TCM cells also acquire de novo CTCF (CCCTC-binding factor) binding sites, which are not only strongly associated with TCM-specific hubs but also linked to increased activities of local gene promoters and enhancers. Specific ablation of CTCF in TCM cells impairs rapid induction of genes in cytotoxic program, energy supplies, transcription, and translation by recall stimulation. Therefore, acquisition of CTCF binding and ChrInt hubs by TCM cells serves as a chromatin architectural basis for their transcriptomic dynamics in primary response and for imprinting the code of "recall readiness" against secondary challenge.

Long non-coding RNA MIR22HG suppresses the chondrogenic differentiation of human adipose-derived stem cells by interacting with CTCF to upregulate CRLF1.

Improved chondrogenic differentiation of mesenchymal stem cells (MSCs) by genetic regulation is a potential method for regenerating articular cartilage. LncRNA MIR22HG has been proven to accelerate osteogenic differentiation, but the regulation mechanism of chondrogenic differentiation is still unclear. Human adipose-derived stem cells (hADSCs) have been widely utilised for bone tissue engineering applications. The present study aimed to examine the effect of MIR22HG on the chondrogenic differentiation of hADSCs. The results confirmed that MIR22HG was downregulated in the process of chondrogenic differentiation. Subsequently, gain- and loss-of-function of MIR22HG experiments showed that the overexpression of MIR22HG suppressed the deposition of cartilage matrix proteoglycans and decreased the expression of cartilage-related markers (e.g. Sox9, ACAN and Col2A1), whereas the knockdown of MIR22HG had the opposite effect. MIR22HG could bind to CTCF (CCCTC-binding factor), and CTCF could bind to the CRLF1 (cytokine receptor-like factor 1) promoter and upregulate CRLF1 gene expression. Besides, inhibition of CRLF1 can reverse the effect of MIR22HG on cell chondrogenic differentiation of hADSCs. Taken together, our outcomes reveal that MIR22HG suppressed chondrogenic differentiation by interaction with CTCF to stabilise CRLF1.

XPF interacts with TOP2B for R-loop processing and DNA looping on actively transcribed genes.

Co-transcriptional RNA-DNA hybrids can not only cause DNA damage threatening genome integrity but also regulate gene activity in a mechanism that remains unclear. Here, we show that the nucleotide excision repair factor XPF interacts with the insulator binding protein CTCF and the cohesin subunits SMC1A and SMC3, leading to R-loop-dependent DNA looping upon transcription activation. To facilitate R-loop processing, XPF interacts and recruits with TOP2B on active gene promoters, leading to double-strand break accumulation and the activation of a DNA damage response. Abrogation of TOP2B leads to the diminished recruitment of XPF, CTCF, and the cohesin subunits to promoters of actively transcribed genes and R-loops and the concurrent impairment of CTCF-mediated DNA looping. Together, our findings disclose an essential role for XPF with TOP2B and the CTCF/cohesin complex in R-loop processing for transcription activation with important ramifications for DNA repair-deficient syndromes associated with transcription-associated DNA damage.

Involvement of CCCTC-binding factor in epigenetic regulation of cancer.

A major global health burden continues to be borne by the complex and multifaceted disease of cancer. Epigenetic changes, which are essential for the emergence and spread of cancer, have drawn a huge amount of attention recently. The CCCTC-binding factor (CTCF), which takes part in a wide range of cellular processes including genomic imprinting, X chromosome inactivation, 3D chromatin architecture, local modifications of histone, and RNA polymerase II-mediated gene transcription, stands out among the diverse array of epigenetic regulators. CTCF not only functions as an architectural protein but also modulates DNA methylation and histone modifications. Epigenetic regulation of cancer has already been the focus of plenty of studies. Understanding the role of CTCF in the cancer epigenetic landscape may lead to the development of novel targeted therapeutic strategies for cancer. CTCF has already earned its status as a tumor suppressor gene by acting like a homeostatic regulator of genome integrity and function. Moreover, CTCF has a direct effect on many important transcriptional regulators that control the cell cycle, apoptosis, senescence, and differentiation. As we learn more about CTCF-mediated epigenetic modifications and transcriptional regulations, the possibility of utilizing CTCF as a diagnostic marker and therapeutic target for cancer will also increase. Thus, the current review intends to promote personalized and precision-based therapeutics for cancer patients by shedding light on the complex interplay between CTCF and epigenetic processes.

CTCF coordinates cell fate specification via orchestrating regulatory hubs with pioneer transcription factors.

CCCTC-binding factor (CTCF), a ubiquitously expressed architectural protein, has emerged as a key regulator of cell identity gene transcription. However, the precise molecular mechanism underlying specialized functions of CTCF remains elusive. Here, we investigate the mechanism through integrative analyses of primary hepatocytes, myocytes, and B cells from mouse and human. We demonstrate that CTCF cooperates with lineage-specific pioneer transcription factors (TFs), including MyoD, FOXA, and PU.1, to control cell identity at 1D and 3D levels. At the 1D level, pioneer TFs facilitate lineage-specific CTCF occupancy via opening chromatin. At the 3D level, CTCF and pioneer TFs form regulatory hubs to govern the expression of cell identity genes. This mechanism is validated using MyoD-null mice, CTCF knockout mice, and CRISPR editing during myogenic differentiation. Collectively, these findings uncover a general mechanism whereby CTCF acts as a cell identity cofactor to control cell identity genes via orchestrating regulatory hubs with pioneer TFs.

Abdominal surgery induces long-lasting changes in expression and binding of CTCF with impact on Major Histocompatibility Complex II transcription in circulating human monocytes.

BACKGROUND: Postoperative immunosuppression has been recognized as an important driver of surgery-related morbidity and mortality. It is characterized by lymphocyte depression and impaired monocyte capability to present foreign antigens to T-cells via Major Histocompatibility Complex, Class II (MHC-II) molecules. In patients with postoperative abdominal sepsis, we previously detected a persisting differential binding of the CCCTC-Binding Factor (CTCF), a superordinate regulator of transcription, inside the MHC-II region with specific impact on human leucocyte antigen (HLA) gene expression. In this prospective exploratory study, we investigated to which extent major surgery affects the MHC-II region of circulating CD14+-monocytes. RESULTS: In non-immunocompromised patients undergoing elective major abdominal surgery, a postoperative loss of monocyte HLA-DR surface receptor density was accompanied by a decline in the transcription levels of the classical MHC-II genes HLA-DRA, HLA-DRB1, HLA-DPA1 and HLA-DPB1. The surgical event decreased the expression of the transcriptional MHC-II regulators CIITA and CTCF and led to a lower CTCF enrichment at an intergenic sequence within the HLA-DR subregion. During the observation period, we found a slow and only incomplete restoration of monocyte HLA-DR surface receptor density as well as a partial recovery of CIITA, HLA-DRA and HLA-DRB1 expression. In contrast, transcription of HLA-DPA1, HLA-DPB1, CTCF and binding of CTCF within the MHC-II remained altered. CONCLUSION: In circulating monocytes, major surgery does not globally affect MHC-II transcription but rather induces specific changes in the expression of selected HLA genes, followed by differential recovery patterns and accompanied by a prolonged reduction of CTCF expression and binding within the MHC-II region. Our results hint toward a long-lasting impact of a major surgical intervention on monocyte functionality, possibly mediated by epigenetic changes that endure the life span of the individual cell.

Corticosteroid-induced chromatin loop dynamics at the FKBP5 gene.

FKBP5 is a 115-kb-long glucocorticoid-inducible gene implicated in psychiatric disorders. To investigate the complexities of chromatin interaction frequencies at the FKBP5 topologically associated domain (TAD), we deployed 15 one-to-all chromatin capture viewpoints near gene promoters, enhancers, introns, and CTCF-loop anchors. This revealed a "one-TAD-one-gene" structure encompassing the FKBP5 promoter and its enhancers. The FKBP5 promoter and its two glucocorticoid-stimulated enhancers roam the entire TAD while displaying subtle cell type-specific interactomes. The FKBP5 TAD consists of two nested CTCF loops that are coordinated by one CTCF site in the eighth intron of FKBP5 and another beyond its polyadenylation site, 61 kb further. Loop extension correlates with transcription increases through the intronic CTCF site. This is efficiently compensated for, since the short loop is restored even under high transcription regimes. The boundaries of the FKBP5 TAD consist of divergent CTCF site patterns, harbor multiple smaller genes, and are resilient to glucocorticoid stimulation. Interestingly, both FKBP5 TAD boundaries harbor H3K27me3-marked heterochromatin blocks that may reinforce them. We propose that cis-acting genetic and epigenetic polymorphisms underlying FKBP5 expression variation are likely to reside within a 240-kb region that consists of the FKBP5 TAD, its left sub-TAD, and both its boundaries.